Tomographic RF Spectroscopy of a Trapped Fermi Gas at Unitarity 



o 
o 

(N 



in 



5-H 

CD 



I 

O 

o 



> ■ 

oo : 

rn ; 

o . 
i> . 
o ■ 



X 



Y. ShinE C. H. Schunck, A. Schirotzek, and W. Ketterle 
Department of Physics, MIT-Harvard Center for Ultracold Atoms, and Research Laboratory of Electronics, 
Massachusetts Institute of Technology, Cambridge, Massachusetts, 02139 
(Dated: February 1, 2008) 

We present spatially resolved radio- frequency spectroscopy of a trapped Fermi gas with resonant 
interactions and observe a spectral gap at low temperatures. The spatial distribution of the spectral 
response of the trapped gas is obtained using in situ phase-contrast imaging and 3D image recon- 
struction. At the lowest temperature, the homogeneous rf spectrum shows an asymmetric excitation 
line shape with a peak at 0.48(4)e_F with respect to the free atomic line, where £f is the local Fermi 
energy. 

PACS numbers: 03.75.Ss, 03.75.Hh, 32.30.Bv 



Ultracold Fermi gases near a Feshbach resonance pro- 
vide new insight into fermionic superfluidity, allowing the 
study of the crossover from Bardeen-Cooper-Schrieffer 
(BCS) superfluids of Cooper pairs to Bose-Einstein con- 
densates (BECs) of dimer molecules. Many aspects of 
the BEC-BCS crossover, including superfluidity have 
been experimentally investigated The properties of 
Fermi gases on resonance at unitarity, where the ampli- 
tude of the scattering length between the fermion atoms 
diverges to infinity and the system shows universal be- 
havior are of great importance to understand the 
crossover physics. Precision measurements of the criti- 
cal temperature S4;] , the interaction energy and collec- 
tive excitations 6] have presented stringent quantitative 
test to the theoretical description of strongly interacting 
Fermi gases. 

Radio-frequency (rf) spectroscopy measures an exci- 
tation spectrum by inducing transitions to different hy- 
perflne spin states. This method has been employed in 
strongly interacting Fermi gases, leading to the obser- 
vation of unitarity limited interactions [3, S|, molecule 
formation on the BEC side of the Feshbach resonance 
as well as pairing in the crossover regime [13, Hi- Rf 
spectroscopy provides valuable information on the na- 
ture of the pairs. Since an rf photon can dissociate bound 
molecules or fermion pairs into the free atom continuum, 
the binding energy of the pairs or the excitation gap is 
determined. Furthermore the excitation line shape is re- 
lated to the wave function of the pairs, e.g. larger pairs 
have narrower lines. However, currently all experimental 
measurements on the excitation spectrum in strongly in- 
teracting Fermi gases 13, 12 1 have been performed with 
samples confined in a harmonic trapping potential so that 
the spectral line shape is broadened due to the inho- 
mogeneous density distribution of the trapped samples, 
preventing a more stringent comparison with theoretical 
predictions 0, 0, E [li ■ 

In this Letter, we demonstrate spatially resolved rf 
spectroscopy of a trapped, population-balanced Fermi 
gas at unitarity at very low temperature. The spatial dis- 
tribution of the rf-induced excited region in the trapped 



gas was recorded with in situ phase-contrast imaging 17 1 
and the local rf spectra were compiled after 3D image re- 
construction. In contrast to the inhomogcncous rf spec- 
trum, the homogeneous local rf spectrum shows a clear 
spectral gap with an asymmetric line shape. We observe 
that the spectral peak shifts by 0.48(4)£i? to higher en- 
ergy and that the spectral gap is 0.30(8)ei:- with respect 
to the free atomic reference line, where ej? is the local 
Fermi energy. This new spectroscopic method overcomes 
the line broadening problem for inhomogeneous samples 
and provides homogeneous rf spectra of a resonantly in- 
teracting Fermi gas revealing the microscopic physics of 
fermion pairs. 

We prepared a degenerate Fermi gas of spin-polarized 
^Li atoms in an optical trap, using laser cooling and 
sympathetic cooling with ^^Na atoms, as described in 
Ref. [l3|- An equal mixture of the two lowest hyperfine 
states |1) and |2) (corresponding to the \F = 1/2, mp = 
1/2) and \F — 1/2, mp = —1/2} states at low magnetic 
field) was created at a magnetic field B — 885 G. A broad 
Feshbach resonance located at Bq — 834 G strongly en- 
hanced the interactions between the two states. The fi- 
nal evaporative cooling by lowering the trap depth and 
all spectroscopic measurements were performed at i? = 
833 G. The total atom number was Nt = 1.0 ± 0.1 x 10^ 
and the radial (axial) trap frequency was fr = 129 Hz 
[fz ~ 23 Hz). The Fermi energy (temperature) of a 
noninteracting equal mixture with the same total atom 
number is Ep = Kf^ f^f/^imtf/^ ^ h x 22.3 kHz 
{Tp = Ep/ks = 1.07 /iK), where h is Planck's constant 
and fcs is Boltzmann's constant. The ratio of the sample 
temperature T to Tp of ~ 0.06 was determined by fitting 
a finite temperature Thomas-Fermi (TF) distribution to 
the whole cloud after expansion. 

Rf spectroscopy was performed by driving atoms in 
state |2} to the next lowest hyperfine state |3) (corre- 
sponding to \F = 3/2, mp — —3/2) at low field) that 
was initially empty. After applying an rf pulse of 1 ms, 
we directly measured the in situ distribution of the den- 
sity difference Ud = fii — n2 in the excited sample, using 
a phase-contrast imaging technique 17|, where rii and 
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FIG. 1: (color online) Radio-frequency (rf) spectroscopy of a 
Fermi gas with in situ phase-contrast imaging, (a) After ap- 
plying an rf pulse, the spatial distribution of the density dif- 
ference between state |1) and |2) is recorded with the phase- 
contrast imaging technique [13]. The density depletion re- 
flects the spin excitation induced by the rf pulse. The dashed 
line indicates the size of the trapped sample, (b) An absorp- 
tion image of an equal mixture without applying an rf pulse. 
The field of view for each image is 205 fim x 680 /im. (c) 
Overall rf spectrum of the inhomogeneous sample is obtained 
by integrating the signal in the phase-contrast images. The 
red line is a Gaussian fit to the spectrum. 



n2 are the densities of atoms in the states |1) and |2), 
respectively. The frequency of the imaging beam was set 
to make the net phase shift in the sample proportional to 
the density difference nd, which was done by zeroing the 
optical signal in an equal mixture without applying an 
rf pulse [17| . Since the initial atom densities in state |1) 
and 1 2) are equal, the density difference nd represents the 
atom number depletion in state |2), the spectral response 

/ M- 

The total spectral response, obtained by integrating 
over the phase contrast images, reproduces earlier re- 
sults 



lol . 11 1 . The phase-contrast images now reveal the 
nature of the observed line shape (Fig. [1]). The spec- 
tral response strongly depends on position. The inner 
region of the cloud, which is at higher density, shows a 
higher resonance frequency. The integrated inhomoge- 
neous spectrum peaks at the rf offset Avi w 10 kHz 20f|. 
The spatially resolved images reveal that at this fre- 
quency, no excitations occur in the center of the cloud, 
but rather in a spatial shell. The rf offset Aiy is measured 
with respect to the resonance frequency of the |2) — |3) 
transition in the absence of atoms in state |1) 21 1. 

Local rf spectra /(r, Ai>) are compiled from the recon- 
structed 3D radial profiles of the density difference. A 
phase-contrast image contains the 2D distribution of the 
column density difference integrated along the imaging 
line, fidix^z) = J nd{r) dy. The excited regions have an 



elliptical shape with the same aspect ratio as the trap, 
A = fz/ fr, confirming the validity of the local density 
approximation. Therefore, we can use elliptically aver- 
aged profiles of the column density difference, nd{r), to 
improve the signal-to-noise ratio, where the ellipse for av- 
eraging is defined as x^ + A^z^ = r^. The 3D radial profile 
Udir) is calculated using the inverse Abel transformation 
of nd{r) [g^l and /(r, Av) cx Udir; Az/). 

With this technique, we obtain homogeneous rf spectra 
as a function of the 3D position, shown in Fig. [21 These 
spectra are the main result of this paper and we now dis- 
cuss their features and implications for our system. The 
local homogeneous rf spectra shows a spectral gap. The 
spectral peak is shifted away from the atomic reference 
line by much more than its line width. Such a gap is not 
observed in the inhomogeneous rf spectrum (Fig. [l^c)) 
where the Gaussian wings overlap with the position of 
the free atomic line. Furthermore, the local rf spec- 
trum reveals an asymmetric line shape of the excitation 
spectrum. For the central region, the peak is located at 
Avp « 15 kHz and the spectral gap, defined as the mini- 
mum energy offset for excitation, is hAvg « /i x 10 kHz. 

The spectral peak position Avp in the local rf spec- 
tra shows a parabolic dependence on the radial position 
(Fig. [Ha)). This can be explained by unitarity, which de- 
mands that all energetic quantities scale with the Fermi 
energy. At unitarity, the only relevant energy scale in the 
system is the Fermi energy ep = h'^i&TT^n )2/3/2to|,[23, 
where n is the atom density in one spin state and m 
is the atomic mass, so the energetic quantities such as 
chemical potential /x and pairing gap energy A are pro- 
portional to i.e., /I = and A = -qep with the 
universal parameters ^ and 77. Therefore the excitation 
spectrum should also scale with the Fermi energy. In an 
external harmonic potential V{r) oc r^, the local Fermi 

energy eF{r)^^i{r)|i={^l^-V{r))/i^eFo{l-r■' I R^), 
where /ig is the global chemical potential, Efo is the local 
Fermi energy at the center, R is the radius of the trapped 
sample and Efq — fJ-o/^ = The spectral peak 

position Ai>p{r) reflects the parabolic radial dependence 
of the local Fermi energy sp^r). 

The local Fermi energy at the center is determined 
from e^o = {R/ Rtp)~'^Ep, where Rtf is the radial 
Thomas- Fermi radius for a noninteracting Fermi gas with 
the same atom number. We obtain Rtp — 67.3 ± 1.1 /xm 
for the measured total atom number and trap frequen- 
cies. The radius of the trapped sample was measured to 
be i? = 56.6 ±1.8 /im, using absorption images like in 
Fig. mb) and fitting the non-saturated outer wing pro- 
file to a zero-temperature TF distribution, giving e_Fo = 
h X 31.5±2.5 kHz. We estimate the universal parameter 

= {R/RtpY = 0.50 ±0.07, which is in goo d agreement 
with previous measurements and 
Quantum Monte Carlo calculations 
in some references). 

The local spectrum at the center of our sample shows 
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FIG. 2: (color online) Spatially resolved rf spectroscopy of a 
trapped Fermi gas. (a) The spectral intensity I{r, Au) was 
obtained from the reconstructed 3D profiles of the density 
difference. See text for the description of the reconstruction 
method. Local rf spectra are shown for (b) r = fim, (c) 
r = 24 iim, and (d) r = 40 fim, whose positions are marked 
by vertical dashed lines in (a). Each spectrum is obtained by 
spatially averaging over 2.5 fim. Rtf is the radial Thomas- 
Fermi radius for a noninteracting Fermi gas with the same 
atom number. The spectral peak position Aj/p(r) in the local 
rf spectra is marked by the black line in (a). The determi- 
nation of Avp is limited to r < 48 /im due to the signal-to- 
noise ratio. The yellow line is a parabolic fit to Avp{r). The 
radius determined by extrapolating the fit to zero rf offset 
is Rp — 53.6 nm indicated by the white down arrow. The 
black up arrow indicates the radius of the trapped sample, 
R = 56.6 /im, measured independently from absorption im- 
ages like Fig. mb). 



the spectral peak at hAi^p = 0.48(4)£i?o and the spec- 
tral gap of hAh'g = 0.30(8)£fo- We determine a radius 
Rp such that Avp{Rp) = 0, extrapolating the parabolic 
curve fit of Avp to zero rf offset (Fig[5][a)). Rp = 53.6 /im 
is slightly smaller than the measured radius R, which we 
attribute to finite temperature effects. Previous studies 
of rf spectroscopy of Fermi gases [13, EH demonstrated 
that the spectral peak shifts to higher energy at lower 
temperature, which is interpreted as the increase of the 
pairing gap energy. In the outer region of lower density, 
the local T/Tp becomes higher, consequently reducing 
hAvp/ep- The observation of Rp being close to R implies 
that our experiment is very close to the zero tempera- 
ture limit. From the relation T/Tpir) oc (1 - r^/i?^)"^ 
we can estimate T/TpiRp) « 15 x T/Tp[Q). If we as- 
sume that the pairing gap energy starts emerging at 
T/Tp w 0.6 [ill, we might infer the local T/Tp < 0.05 at 
the center, close to our measured temperature. Although 
hAvp/sp is almost constant over the whole sample, the 
line width increases in the outer region. 

The homogeneous rf spectra measured in our experi- 
ment allow a direct comparison with theoretical predic- 
tions. However, a comprehensive theoretical interpreta- 
tion of the rf spectrum of a strongly interacting Fermi gas 
is not available yet and we discuss here an extrapolation 
of BCS theory to strong unitarity limited interactions. 
Rf spectroscopy measures a single-particle spin excita- 
tion spectrum, since an rf photon changes the spin state 
while imparting negligible momentum. The conventional 
picture of rf spectroscopy of pairs is a photodissociation 
process: the initial |1) — |2) bound state, which can be 
molecules or fcrmion pairs, breaks into free particles in 
state |1) and |3). In a BCS superfluid, the free particle 
in state |1) is regarded as a quasiparticle, so after the 
spin transition the whole system can be described as the 
excited BCS state with one quasiparticle and one free 
particle in state |3). With the assumption of no inter- 
actions between state |1) and |3), the rf photon energy 
offset would be hAiy — E-k — M + Ji^k"^ /2m, where the 

first term E^k ~ \J A"^ + {h^k'^ /2m — iiY is the energy 
cost for generating one quasiparticle excitation with mo- 
mentum —k, the second term is for removing one atom 
in state |2), and the last term is the kinetic energy of the 
atom in state |3) with momentum k. 

The measured FWHM line width is about two times 
narrower than predicted by the simple model described 
above. The model spectrum shows a very long tail cor- 
responding to high momentum contributions. This dis- 
crepancy might be due to modification of the BCS ex- 
pressions in the unitarity regime. The narrow peak might 
imply that the pair wave function is narrower in momen- 
tum space and therefore more spatially extended than 
the BCS prediction. Another extension of our simple 
model should address the interactions between atoms in 
state |1) and |3). The mean-field interaction energies due 
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to |1) — |2) and |1) — |3) interactions have been empir- 
ically assumed to have the same unitarity limited value 
because of the proximity of a |1) — |3) Feshbach reso- 
nance at B=690 G 1^, T^. The recent experiments with 
imbalanced mixtures and higher densities showed some 
deviations from this assumption flT|. 

A localized spin excitation, induced by an rf pulse, 
eventually diffuses over the sample. This ultimately lim- 
its the pulse duration and therefore the spatial resolution. 
Using phase-contrast imaging, we monitored the evolu- 
tion of the spatial structure in the excited sample with 
various delay time after applying an rf pulse. The shell 
structure was well preserved even after 5 ms and only 
some broadening was observed, showing that during the 
1-nis pulse the dynamic evolution of the density differ- 
ence profiles is negligible. 

We found that atoms in state |3) rapidly decayed in 
the sample, although the total atom number difference 
between state |1) and |2) did not change over time. At 
the center of the sample where the total atom density is 
about 8 X 10^^ cm~^, the life time of atoms in state |3) was 
measured to be less than 0.4 ms. Since a |1) — |3) mixture 
and a |2) — |3) mixture are stable at B = 833 G [33 |. 
the decay should be associated with exoergic molecule 
formation via three-body collisions involving one atom 
from each spin state. We observed that the loss of one 
atom in state |3) was accompanied by loss of one atom 
in state |1) and one atom in state |2), supporting the 
three-body loss mechanism 33 1. 

Our new technique of spatially resolved rf spectroscopy 
should be able to address the important questions also at 
finite temperature. One question is whether the observed 
double peak structure [iQ, [lH of an atomic line and a 
pairing peak is purely inhomogeneous, or whether it is 
possible to have local coexistence of paired and unpaired 
atoms. This is of course possible on the BEG side of 
the Feshbach resonance where in a certain temperature 
range, bound molecules and thermally dissociated free 
atoms coexist, but it is an open question, how this picture 
will change in the BEC-BGS crossover. 

In conclusion, we present spatially resolved rf spec- 
troscopy of a trapped Fermi gas, using the in situ phase- 
contrast imaging technique. The homogeneous rf spec- 
tra of a Fermi gas at unitarity provide a benchmark for 
a complete theoretical description, which should reveal 
microscopic details of the paired states. 
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manuscript. This work was supported by the NSF and 
ONR. 
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